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Computational  Fluid  Dynamics  Study  for  a  Deep  Stall  Air  Vehicle 
INTRODUCTION 

Unmanned  air  vehicles  (UAVs)  have  become  an  integral  part  of  the  battle-space 
environment.  Small  camera-equipped  vehicles  have  shown  their  value  in  recent  conflicts.  These 
vehicles  can  be  assembled  and  launched  from  small  spaces,  but  recovery  of  these  vehicles  still 
poses  difficulties.  Recovery  of  these  aircraft  typically  requires  a  long  open  space  for  landing.  In 
order  to  minimize  the  space  required  for  recovery  it  is  desired  to  design  a  vehicle  capable  of 
flight  at  high  angles  of  attack.  High  angle-of-attack  flight  could  also  provide  greater  obstacle 
avoidance  capabilities  when  flying  in  urban  terrain  since  the  turning  radius  in  high  angle-of- 
attack  flight  is  greatly  reduced  as  compared  to  that  of  conventional  forward  flight.  Sensing 
capabilities  could  also  be  enhanced  through  the  use  of  sensor  pointing  without  the  complexity 
and  weight  of  morphing  aircraft  or  gimbaled  sensors. 

Small  UAVs  have  received  recent  attention  for  several  urban  missions  due  to  their  rapidly 
maturing  flight  capabilities.  These  small  man-portable  vehicles  feature  a  deep-stall  mode  in 
which  the  elevator  is  deflected  and  the  power  is  reduced.  The  aircracft  then  descends  steeply  in  a 
near  level  attitude  until  it  touches  the  ground.  An  example  of  such  a  vehicle  is  the  Raven, 
developed  by  Aerovironment  [1],  which  utilizes  the  wing  geometry  to  remain  stable  in  deep- 
stall.  At  high  angles  of  attack,  the  flight  of  these  vehicles  are  often  characterized  by  wing  rock. 
A  possible  cause  for  this  phenomenon  is  the  unsteady  shedding  of  vortices  from  the  wing  and 
further  interaction  of  this  flow  with  the  airframe  downstream. 

In  order  to  operate  these  small  UAVs  in  deep-stall,  Aaron  Kahn  of  the  U.S.  Naval  Research 
Laboratory  has  developed  an  adaptive  controller  [2]  capable  of  regulating  vehicle  attitude  in  the 
presence  of  large  model  errors.  The  focus  of  the  current  study  is  to  obtain  the  aerodynamic 
characteristics  of  the  vehicle  in  the  deep  stall  regime.  An  unstructured  grid-based  unsteady 
Navier-Stokes  solver  with  automatic  adaptive  remeshing  capability,  called  FEFLO,  was  used  to 
compute  the  flow  about  the  deep-stall  vehicle. 


THE  INCOMPRESSIBLE  FLOW  SOLVER 

The  governing  equations  employed  are  the  incompressible  Navier-Stokes  equations  in 
Arbitrary  Lagrangian-Eulerian  (ALE)  formulation  which  are  written  as 

d\ 

—  +vfl-Vv  +  Vp  =  V-cr,  (1) 

dt 

V*  v  =0,  (2) 

where  p  denotes  the  pressure,  \a  -  v  -  w  the  advective  velocity  vector,  where  v  is  the  flow 
velocity  and  w  is  the  mesh  velocity  w  and  the  material  derivative  is  with  respect  to  the  mesh 
velocity  w.  Both  the  pressure  p  and  the  stress  tensor  o  have  been  normalized  by  the  (constant) 
density  p  and  are  discretized  in  time  using  an  implicit  time  stepping  procedure.  Thus  the 
equations  are  Eulerian  for  zero  mesh  velocity  and  Lagrangian  if  the  mesh  velocity  is  the  same  as 
the  flow  velocity.  The  present  time-accurate  flow  solver  is  discretized  in  space  using  a  Galerkin 
procedure  with  linear  tetrahedral  elements.  The  details  of  the  flow  solver  have  already  been 
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discussed  extensively  elsewhere,  Ramamurti  et.  al.  [3-6],  in  connection  with  successfully 
validated  solutions  for  numerous  2-D  and  3-D,  laminar  and  turbulent,  steady  and  unsteady  flow 
problems. 

THE  COMPUTATIONAL  MODEL 

The  initial  critical  and  time  consuming  step  for  flow  simulations  is  the  generation  of  a  3-D 
surface  model  that  is  suitable  for  computational  fluid  dynamics.  This  is  achieved  by  converting 
the  longitudinal  cross-sectional  drawing  of  the  vehicle  to  a  digitized  format  using  ImageJ  [7], 
For  example,  Fig.  1  shows  the  drawing  of  the  wing  tip  and  the  digitized  points.  These  digitized 
sections  are  then  assembled  using  FECAD,  the  CAD  software  used  in  conjunction  with  the 
FEFLO  software  suite.  The  CAD  representation  of  the  complete  deep-stall  vehicle  is  shown  in 
Fig.  2. 


Fig.  1.  Digitization  of  the  vehicle  geometry:  a)  drawing  of  the  wing  tip  and  b)  digitized  points. 


Fig.  2.  CAD  representation  of  the  Deep-Stall  Vehicle. 


Computational  Fluid  Dynamics  Study  for  a  Deep  Stall  Air  Vehicle 


3 


RESULTS  AND  DISCUSSION 

In  order  to  determine  the  steady  state  aerodynamic  characteristics  of  this  vehicle,  a 
tetrahedral  mesh  consisting  of  1.2  M  points  and  6.8  M  tetrahedral  elements.  The  geometry  used 
to  evaluate  the  pitch  characteristics  employed  a  symmetry  plane  and  therefore  only  half  the 
geometry  shown  in  Fig.  2  is  used.  For  determining  side-slip  characteristics  and  computing 
unsteady  flow  past  the  vehicle,  a  domain  consisting  of  approximately,  2.8  M  points  and  16.0  M 
tetrahedral  elements  is  employed.  The  forces  and  moments  were  obtained  by  integrating  the 
surface  pressure  over  the  various  segments  of  the  vehicle,  such  as  the  wing  and  the  fuselage,  in 
the  fixed  coordinate  system  (XYZ),  shown  in  Fig.  3.  The  resulting  forces  were  transformed  to 
the  vehicle  coordinate  system  (xyz)  using  the  rotation  transformation,  Eq.  1.  In  this  coordinate 
system,  the  x-axis  is  aligned  with  the  flow  velocity  vector,  and  a  and  (3  are  the  angle  of  attack 
and  side-slip,  respectively.  In  order  to  keep  the  transformation  general,  the  roll  angle  (J)  is  also 
included. 


Y 


Fig.  3.  Computational  Coordinate  system  for  Deep-Stall  Vehicle. 
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The  computed  forces  (Fx,  FY  ,FZ)  are  transformed  to  their  respective  drag,  lift  and  side  forces 
(D,  L,  N)  which  are  normalized  using  the  free  stream  velocity  V,  the  density  of  air  p,  and  the  area 
of  the  wing  Aw. 
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The  moments  are  further  normalized  by  the  average  wingspan,  S  =  AW/C,  where  C  =  0.171m  is 
the  chord  length  of  the  wing. 

Effect  of  angle  of  attack: 

The  steady  flow  characteristics  are  obtained  for  a  vehicle  speed,  V  =  13.41  m/s,  and  for 
various  angles  of  attack  and  side-slip.  Figure  4  shows  the  lift,  drag  and  pitch  moment  variation 
as  the  angle  of  attack  is  increased.  The  forces  and  moments  on  just  the  wing  surfaces  and  the 
complete  vehicle  are  shown.  It  can  be  seen  that  for  a’s  in  the  range  of  (-10°,  15°),  the  lift 
variation  is  nearly  linear,  as  expected.  The  vehicle  stalls  for  a  >  15°  shown  by  the  drop  in  lift 
and  shows  a  modest  increase  in  lift  beyond  25°.  At  a  =  0°,  the  drag  force  on  the  vehicle  is 
approximately  0.29  N,  yielding  CD0  of  4.614  x  10"3;  the  corresponding  lift  is  8.48  N  and  CL  is 
0.2392.  At  a  =  15°,  the  lift  and  drag  forces  are  58.878  N  and  5.2  N,  respectively.  The  Reynolds 
number  for  this  flow  based  on  the  chord  length  is  approximately  153,000.  Hence,  a  viscous  flow 
past  this  configuration  was  simulated  at  both  0°  and  15°  angles  of  attack.  The  resulting  forces 
show  that  at  a  =  0°,  the  lift  and  drag  forces  are  8.29N  and  0.62N,  respectively;  at  a  =  15°,  the  lift 
and  drag  forces  are  57.97  N  and  5.55  N,  respectively. 


b.  Drag  Coefficient  c.  Pitch  moment  Coefficient 

Fig.  4.  Lift,  drag  and  pitch  moment  characteristics  of  the  Deep-Stall  Vehicle. 
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Effect  of  Side  Slip  angle: 

In  order  to  assess  the  characteristics  of  the  vehicle  with  respect  to  the  side  slip  angle,  the  flow 
past  the  vehicle  was  simulated  over  a  range  ±15°,  with  the  vehicle  trimmed  at  two  different 
angles,  0°  and  3°.  Figures  5-8  show  the  variation  of  lift,  drag  and  side  forces  and  the  moment 
coefficients.  These  results  show  that  the  side  force  coefficient,  CN,  and  the  moment  coefficients 
about  the  x  and  y-axes  are  linear  with  respect  to  the  side-slip  angle,  as  expected.  The  lift,  drag 
and  roll  moment  coefficients  show  a  parabolic  variation.  As  the  trim  angle  is  increased  form  0° 
to  3°,  the  CL  increases  by  0.3,  Figs.  5a  and  7a,  and  the  CD  also  increases,  while  the  CN  remains 
nearly  unchanged,  as  shown  in  Figs.  5c  and  7c. 


-1.0 


20.0 


a.  Lift  b.  Drag  c.  Side  force 

Fig.  5.  Lift,  drag  and  side  force  coefficients,  vehicle  trim  angle  =  0°. 


Fig.  6.  Moment  Coefficients,  vehicle  trim  angle  =  0°. 


P  (deg.)  p  (deg.) 


P  (deg.) 


a.  Lift  b.  Drag  c.  Side  force 

Fig.  7.  Lift,  drag  and  side  force  coefficients,  vehicle  trim  angle  =  3°. 
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P  (deg.)  p  (deg.)  p  (deg.) 


a.  C 


Mx 


b.  CMy 

Fig.  8.  Moment  Coefficients,  vehicle  trim  angle  =  3°. 
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Effect  of  Rudder  deflection: 

The  effect  of  the  rudder  deflection  on  the  forces  and  moments  were  computed  for  two 
deflections  of  ±5°,  and  for  three  side  slip  angles  of  0°  and  ±15°.  Figures  9a  and  b  show  that  CL 
variation  with  respect  to  the  rudder  deflection  remains  same  while  the  CD  for  rudder  deflection 
6  =  5°  and  |3  =  -15°  is  the  same  as  that  of  6  =  -5°  and  (3  =  15°,  as  expected.  The  side  force 
coefficient  is  increased  by  nearly  0.023,  as  shown  in  Fig.  9c.  The  change  in  the  yaw  moment, 
CMy,  due  to  the  rudder  deflection  of  10°  is  nearly  0.012,  and  the  change  in  the  roll  moment,  CMx, 
0.26 


0.24 

& 

0.22 


0.20  _ 

-20.0  -10.0  0.0  10.0  20.0  -20.0  -10.0  0.0  10.0  20.0  -20.0  -10.0  0.0  10.0  20.0 

P  (deg.)  p  (deg.)  p  (deg.) 


a.  Lift  b.  Drag  c.  Side  force 

Fig.  9.  Lift,  drag  and  side  force  coefficients,  vehicle  trim  angle  =  0°. 


Fig.  10.  Moment  Coefficients,  vehicle  trim  angle  =  0°. 
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is  0.001,  as  shown  in  Figs.  10a  and  b.  The  pitch  moment  remains  unchanged  due  to  the  rudder 
deflection,  as  expected.  Similar  results  are  obtained  for  the  3°  trim  case,  in  Figs.  11  and  12. 


P  (deg.)  p  (deg.) 


-20.0  -10.0  0.0  10.0  20.0 

P  (deg.) 


a.  Lift  b.  Drag  c.  Side  force 

Fig.  11.  Lift,  drag  and  side  force  coefficients,  vehicle  trim  angle  =  3°. 


Fig.  12.  Moment  Coefficients,  vehicle  trim  angle  =  3°. 


Effect  of  Aileron  deflection: 

In  order  or  study  the  effect  of  the  aileron  deflection  on  the  roll  characteristics  of  the  vehicle, 
the  left  aileron  is  deflected  in  a  counterclockwise  (+ve  0)  direction  about  the  Z-axis  (Fig. 2),  and 
the  right  aileron  is  deflected  in  a  clockwise  direction  by  the  same  angle.  Simulations  are 
performed  for  four  aileron  deflections  and  for  side-slip  angles  of  0°  and  15°.  Figure  13a  shows 


0(deg.)  o(deg.) 


8  (deg.) 


a.  Lift  b.  Drag  c.  Side  force 

Fig.  13.  Lift,  drag  and  side  force  coefficients,  vehicle  trim  angle  =  0°. 
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that  the  lift  remains  nearly  unchanged  for  |3  =  0°  and  a  loss  of  lift  for  |3  =  15°.  The  drag  force 
increases  as  the  aileron  deflection  is  increased,  Fig.  13b,  and  the  side  force  remains  constant  for 
non-zero  side-slip  angle,  Fig.  13c.  The  effect  of  the  aileron  deflection  on  the  moments  is  shown 
in  Fig.  14.  The  roll  moment  and  the  pitch  moment  for  non-zero  side  slip  varies  linearly  with  the 
aileron  deflection,  shown  in  Figs.  14a  and  c,  while  the  yaw  moment  and  the  pitch  moment  at  zero 
side  slip  are  nearly  constant,  shown  in  Figs.  14b  and  c.  Similar  results  are  obtained  for  the 
vehicle  trimmed  at  3°,  Figs.  15  and  16. 


Fig.  14.  Moment  Coefficients,  vehicle  trim  angle  =  0°. 


a.  Lift  b.  Drag  c.  Side  force 

Fig.  15.  Lift,  drag  and  side  force  coefficients,  vehicle  trim  angle  =  3°. 
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b.  CMy 

Fig.  16.  Moment  Coefficients,  vehicle  trim  angle  =  3°. 
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Effect  of  the  Elevator: 

The  trailing  32%  of  the  tail  section  is  used  as  the  elevator.  This  section  is  deflected  Se  =  ±5°, 
and  for  the  side-slip  angles  of  ±15°,  for  both  trim  conditions.  The  results  show  that  CL  and  CMz, 
Fig.  17,  exhibit  the  maximum  variation  with  respect  to  the  elevator  deflections,  while  all  the 
other  coefficients  remain  nearly  constant. 
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Fig.  17.  Lift,  drag  and  pitch  moment  coefficients,  vehicle  trim  angle  =  0°. 
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Fig.  18.  Lift,  drag  and  pitch  moment  coefficients,  vehicle  trim  angle  =  3°. 


Deep-Stall  Flow: 

The  vehicle  configuration  was  modified  for  deep-stall  flows  by  rotating  the  tail  about  the  1/3 
chord  location  by  70°.  The  unsteady  flow  past  this  configuration  was  simulated  at  an  angle  of 
attack  a  =  60°  and  70°  and  an  inflow  velocity  V  =  1.67  m/s.  Figures  19a  and  b  show  the 
unsteady  forces  in  the  coordinate  system  xyz,  Fig.  2,  where  x-axis  is  aligned  with  the  inflow 
velocity.  For  a  =  60°,  the  forces  were  averaged  in  the  interval  t  =  4.0s  and  4.5s.  The  mean  drag 
force  on  the  vehicle  is  0.734  N  and  the  mean  lift  force  is  0.412  N.  The  mean  pitch  moment  about 
the  origin  is  0.249  N-m.  For  a  =  70°,  the  forces  and  moments  were  averaged  in  the  interval  t  = 
7.0s  to  7.5s.  The  mean  drag  force  is  0.85  N  and  the  mean  lift  force  is  0.36  N.  The  mean  pitch 
moment  is  0.3  N-m.  Figure  20  shows  the  pressure  distribution  on  the  surface  of  the  vehicle  at 
a  =  70°. 


F  (N) 
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time  (sec.)  time  (sec.) 


a.  a  =  60°  b.  a  =  70° 

Fig.  19.  Unsteady  forces  for  a  deep  stall  configuration,  V  =  1.67  m/s. 


l-io.o 


b.  bottom  view 

Fig.  20.  Surface  pressure  distribution  on  a  Deep  Stall  vehicle,  V  =  1.67  m/s,  a  =  70°. 
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The  deep-stall  characteristics  are  also  obtained  for  a  velocity  V  =  5.79  m/s,  for  a  range  of  a's 
from  50°  to  90°.  The  unsteady  forces  and  moments  were  averaged  over  0.5s  after  the  initial 
transients  have  subsided  and  are  shown  in  Fig.  21a  and  b.  The  drag,  lift  and  pitch-moments  are 
piecewise  linear  with  two  linear  segments  between  50°  -  70°  and  70°  -  90°.  Figure  22a  and  b 
show  the  variation  of  the  mean  forces  and  moments  with  side-slip  angle  is  nearly  linear  for  the 
side  force  an  the  roll  and  yaw  moments. 


Fig.  21.  Mean  forces  and  moments  for  Deep-Stall  vehicle,  V  =  5.79  m/s,  (3  =  0°. 


Fig.  22.  Mean  forces  and  moments  for  Deep-Stall  vehicle,  V  =  5.79  m/s,  a  =  70°. 

Next,  the  roll  characteristics  of  the  vehicle  under  deep-stall  conditions  were  determined  by 
deflecting  the  left  aileron  through  ±20°  deflections.  The  results  show  that  the  mean  forces  and 
moments  remain  nearly  constant  as  shown  in  Fig.  23a  and  b.  The  effect  of  rudder  deflection  on 
the  vehicle  was  studied  for  ±5°  deflections  of  the  rudder.  The  results  show  that  the  fx  =  10.4  N, 
fy  =  4.38  N  and  Mz  =  3.67  N-m.  The  yaw  moment  varies  from  0.027  N-m  to  0.019  N-m  as  the 
rudder  deflection  is  changed  form  -5°  to  +5°. 
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Fig.  23.  Variation  of  Mean  forces  and  moments  for  Deep-Stall  vehicle  with  left  aileron 
deflection,  V  =  5.79  m/s,  a  =  70°. 

All  of  the  above  results  were  used  to  obtain  the  control  derivatives  used  by  Aaron  Kahn  [2] 
for  the  controller  development.  The  unsteady  trajectory  of  the  vehicle  is  next  computed  using 
the  vehicle  properties  such  as  mass,  center  of  gravity  and  the  moments  of  inertia  about  the 
principal  axes  of  the  vehicle.  Initial  simulations  were  performed  at  V  =  5.79  m/s  at  a  =  70°,  with 
the  vehicle  constrained  in  pitch  and  roll.  These  results  coupled  with  the  controller  simulations  of 
Aaron  Kahn  showed  that  the  descent  rate  of  the  vehicle  would  be  too  fast  and  not  feasible. 
Hence,  the  configuration  for  the  deep-stall  recovery  trajectory  simulation  was  modified  by 
setting  the  free  stream  velocity  to  V  =  10.97  m/s  at  a  =  40°.  The  tail  was  rotated  to  align  with 
the  in  flow  angle  and  both  the  ailerons  were  deflected  by  15°.  In  addition,  the  vehicle  was 
initially  pitched  up  by  7.15°.  The  unsteady  flow  past  this  configuration  was  obtained  with  the 
trajectory  of  the  vehicle  constrained  to  pitch  and  roll  due  to  the  computed  forces  and  moments. 
Also,  the  vehicle  was  constrained  in  all  the  translational  degrees  of  freedom.  The  computed 


Fig.  24.  Time  history  of  forces  and  moments  on  a  deep-stall  vehicle,  V  =  10.97  m/s,  a  =  40°. 
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Fig.  25.  Trajectory  of  the  deep-stall  vehicle,  V  =  10.97  m/s,  a  =  40°. 


Fig.  26.  Surface  pressure  distribution  on  the  vehicle  undergoing  deep-stall  maneuver,  t  =  2.38  s. 

unsteady  forces  and  moments  are  shown  in  Fig.  24.  The  pitch  and  roll  motions  of  the  vehicle  is 
shown  in  Fig.  25.  The  results  show  that  the  vehicle  undergoes  initial  pitch  oscillations  and  settles 
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into  a  final  pitch  down  attitude  of  nearly  6°,  while  it  rolls  initially  and  stabilizes  later.  The 
motion  is  consistent  with  the  findings  of  the  reduced  model  simulation  of  Aaron  Kahn  shown  in 
Fig.  25.  The  instantaneous  pressure  distribution  on  the  surface  of  the  vehicle,  Figs.  26a  and  b 
show  the  unsteadiness  on  the  top  surface  of  the  wing  and  near  stagnation  pressure  on  the  bottom 
of  the  wing  and  the  fuselage.  The  unconstrained  trajectory  of  the  pitch  and  roll  motions  of  this 
vehicle  with  differential  aileron  deflections  were  also  computed.  Two  cases  of  aileron 
deflections  were  considered.  The  first  one  was  with  the  right  aileron  deflected  down  by  15° 
while  the  left  aileron  is  flush  with  the  wing.  For  this  case,  the  unsteady  forces  and  moments  are 
shown  in  Fig.  27  and  the  pitch  and  roll  trajectory  is  shown  in  Fig.  28.  The  pitch  motion  is 
similar  in  trend  to  the  previous  case  while  the  roll  is  oscillatory  about  positive  angle  of  nearly  4°, 
as  expected.  The  second  case  was  simulated  with  the  right  aileron  deflected  up  by  30°  while  the 


Fig.  27.  Time  history  of  forces  and  moments  on  a  deep-stall  vehicle  with  differential  aileron 
deflection  A0  =  15°,  V  =  10.97  m/s,  a  =  40°. 


Fig.  28.  Trajectory  of  a  deep-stall  vehicle  with  differential  aileron  deflection  A0  =  15°, 
V  =  10.97  m/s,  a  =  40°. 
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left  aileron  is  flush  with  the  wing.  The  unsteady  force  and  moments  are  shown  in  Fig.  29  and  the 
trajectory  of  the  vehicle  is  shown  in  Fig.  30.  After  4s  of  simulation,  the  vehicle  achieves  a  nearly 
steady  downward  pitch  angle  of  4.75°  and  while  the  roll  has  not  achieved  a  quasi  steady  state. 


Fig.  29.  Time  history  of  forces  and  moments  on  a  deep-stall  vehicle  with  differential  aileron 
deflection  A0  =  30°,  V  =  10.97  m/s,  a  =  40°. 


Fig.  30.  Trajectory  of  a  deep-stall  vehicle  with  differential  aileron  deflection  A0  =  30°, 
V  =  10.97  m/s,  a  =  40°. 
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SUMMARY 

The  3-D  steady  and  unsteady  computations  of  the  flow  past  a  deep-stall  air  vehicle  were 
performed.  The  aerodynamic  characteristics  of  the  vehicle  were  mapped  out  for  a  wide  range  of 
angles  of  attack,  side-slip  and  control  surface  deflections.  The  unsteady  characteristics  of  this 
vehicle  were  also  computed  in  the  deep-stall  regime.  These  results  were  used  for  the 
development  of  dynamic  model  of  the  controller.  For  the  deep-stall  configurations,  the  unsteady 
forces  and  moments  were  computed  by  the  integrating  the  surface  distributions  and  the 
unconstrained  trajectory  of  the  vehicle  is  predicted  based  on  the  6-dof  model. 
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